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Abstract 

It is known that the permittivity of the ferroelectric films is affected by several phenomena, which 
deteriorate the material quality (e.g. the redistribution of the crystal lattice defects, appearance 
of the electrode- adjacent non-ferroelectric layers or the spontaneous polarization screening due 
to a free charge injection across the electrode-adjacent layer, etc.). It is also known that the 
permittivity of ferroelectric polydomain films is controlled by the sum of two contributions: the 
crystal lattice (intrinsic) contribution and the domain wall movement (extrinsic) contribution. It 
is the latter one, which is very sensitive to the aforementioned phenomena and which plays a key 
role in the deterioration of the dielectric response of the ferroelectric polydomain films. In this 
Article, there is presented a method for the identification of the process, which is responsible for 
the ferroelectric ageing. The method is based on the analysis of the evolution of both the linear and 
nonlinear permittivity during ageing. Applicability of the method is theoretically demonstrated 
on four ageing scenarios in two qualitatively different systems where the evolution of the nonlinear 
permittivity is controlled, first, by a redistribution of the pinning centers on the domain wall 
and, second, by microstructural changes at the interface between the ferroelectric layer and the 
electrode. It is shown that each ageing scenario is characterized by unique trend in the evolution 
between the linear and nonlinear part of the permittivity, which can be verified experimentally. 
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I. INTRODUCTION 



Ferroelectric materials play more and more important role in many contemporary elec- 
tronic devices. Due to the fundamental thermodynamic reasons, researchers and applied 
scientists integrating the ferroelectric materials into silicon-based semiconductor chips often 
face problems with a severe deterioration of the quality of the integrated ferroelectric films. 
These deterioration problems may manifest themselves during characterization experiments 
such as ferroelectric hysteresis loop measurements, pulse switching curve measurements, or 
small signal permittivity measurements, and cause phenomena called imprint, fatigue, age- 
ing, etc. Unfortunately, determination of the "microstructural" mechanism predominantly 
responsible for an observed deteriorating effect in a particular ferroelectric sample is always 
very difficult, since several different microstructural mechanisms yield the similar manifesta- 
tions in the aforementioned experiments. In order to get a better insight into microstructural 
mechanisms in ferroelectric samples, which would serve as a feedback for the optimization of 
the fabrication processes, a comprehensive analysis of the results of several characterization 
techniques should be performed. 

The aforementioned issues have motivated the study presented in this work, where we 
analyze the possibility to use the measurements of nonlinear permittivity for the determi- 
nation of mechanisms responsible for the ageing of polydomain ferroelectrics. The principle 
of our study is based on the fact that the dielectric response of the polydomain ferroelectric 
film is controlled by the sum of two contributions: the crystal lattice (intrinsic) contribution 
and the domain wall movement (extrinsic) contribution. Since the latter one is known to 
be sensitive to the geometry of the domain pattern and to several microstructural features 
of the material such as crystal lattice defects or free charges, which may work as pinning 
centers and may reduce the mobility of the domain walls, it is very likely that a compre- 
hensive analysis of nonlinear dielectric permittivity may offer a simple and useful tool for 
getting some additional information about the microstructural processes responsible for the 
ferroelectric ageing. 

The concept of the ferroelectric domain boundary contribution to permittivity has been 
discussed for the first time by Kittel in 1951 [1]. In the following decade, the existence 
of extrinsic contributions has been discussed several times in the literature, mainly in the 
context of Kittel's theory on the resonance character of permittivity of polydomain ferro- 
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electric in the microwave frequency region due to effective domain wall inertia [2, 3]. In 
I960, Misafova [4] made an early discovery of putting the existence of extrinsic contribu- 
tions to permittivity in connection with ferroelectric ageing. Nevertheless, a definite answer 
was not given to the question of the magnitude of the extrinsic contribution with respect to 
the intrinsic one until the work by Fousek [5] where the methodology of the experimental 
determination of the extrinsic and intrinsic contributions has been established in a sense 
understood today, i.e. the intrinsic contribution corresponds to the dielectric response of 
the poled single-domain sample and the extra enhancement in the dielectric response, which 
is measured on a polydomain sample, corresponds to the extrinsic contribution. After this 
breakthrough in understanding of the concept, the extrinsic contributions to permittivity, 
piezoelectric coefficients, and elastic compliance became the subject of an intense theoretical 
and experimental research [6, 7, 8, 9]. Up to date, the extrinsic contributions have been con- 
sidered something like a material curiosity of polydomain samples, which might be, however, 
beneficially used in applications. Nevertheless, the substantial drawback of the application 
of polydomain ferroelectric in electronic devices is the essential instability of the domain 
pattern, which is often manifested by a severe ageing of dielectric properties. 

In this work, we present a method based on the comprehensive analysis of the evolution 
of the nonlinear dielectric response, which is caused by an evolution in the domain pattern 
microstructure. In Sec. II we analyze the intrinsic and extrinsic contributions to the nonlin- 
ear dielectric response of the ferroelectric polydomain film in detail. Section III presents a 
general method to identify the dominant microstructural phenomenon, which is responsible 
for the ageing of the ferroelectric. In Sec. IV we demonstrate the application of our method 
on two qualitatively different systems. First, the system where the extrinsic permittivity is 
controlled by the bending movements of pinned domain walls (Sec. IV A) and, second, the 
system with electrode-adjacent passive layers (Sec. IV B). We show that both systems under 
consideration may manifest several microstructurally different ageing scenarios, which are 
characterized by different features of the evolution of nonlinear permittivity. 
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FIG. 1: Scheme of the ferroelectric capacitor of thickness h with the electric field E applied to the 
material, which is in the polydomain state. Symbols «W and v ^ stand for the volume fractions 
of the ferroelectric, which is in the polarization state and P^ 2 \ respectively. 

II. NONLINEAR DIELECTRIC RESPONSE OF THE POLYDOMAIN FERRO- 
ELECTRIC 

At first we analyze in detail the intrinsic and extrinsic contribution to the nonlinear di- 
electric response of the polydomain ferroelectric film, which is shown in Fig. 1. We consider 
a ferroelectric capacitor of thickness h with the electric field E applied to the material, which 
is split by 180° domain walls into the lamellar pattern of antiparallel domains at the two 
states with polarizations denoted by the symbols P/ 1 ^ and P/ 2 ^ . The volume fractions of 
the domains are denoted by the symbols and v^ 2 \ In order to calculate the polarization 
response of the crystal lattice of the ferroelectric film, we adopt the incomplete thermody- 
namic potential (per unit volume of the sample) g(Pi, Ei) = $(Pi) — P%Ei, which includes, 
first, the energy associated with the polarization of the crystal lattice $(P;), and, second, the 
subtracted work produced by external voltage sources —PiEi. In order to present the basic 
principles of our method, we limit our considerations to the simplest case of the uniaxial 
ferroelectric with the vectors of the electric field Ei and the polarization Pi perpendicular 
to the top and bottom electrodes of the capacitor. A generalization to a more complicated 
system is more or less straightforward. For the simplification of the notation, we will omit 
the indices "i" in the symbols for the vectors Pj and Ei in the following text. Then, the 
three lowest terms in the expansion of the function g(P, E) with respect to the polarization 
P of the crystal lattice equals: 



g(P, E) = %+ X -aP 2 + - A (3P A - PE. 



(1) 
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The polarization response of the crystal lattice can be found from the condition for the 
minimum of the thermodynamic potential g: 



dg(P,E)/dP = 0. (2) 

In the ferroelectric phase, when a < and /3 > 0, the polarization response in the two stable 
polarization states of the crystal lattice in the adjacent domains can be expressed in a form 
of Taylor series with respect to the electric field P: 

P«(P) « P + e oXc E - P* + p3 _ . . . , (3a) 

P^\E) » -P + eoXc£ + ^ ^ 2 + ^ ^ 3 + • ■ ■ , (3b) 



where Po = y— /3 /a is the spontaneous polarization and Xc = — l/(2^o tt ) is the dielectric 
susceptibility of the crystal lattice. 

Now we can consider that the average polarization of the polydomain ferroelectric film P/ 
is given by the volume weighted average of the polarization response in each domain state: 

P f (E) = v w P (1) (P) + v {2) P (2) (P), (4) 

where + = 1. Substituting Eqs. (3) into Eq. (4), one gets: 

P f (E) = («W - v^) (p - ^ P 2 ) + .oXei? + ^ P 3 + • • • . (5) 

It is convenient to express the above expression in terms of the net spontaneous polarization 
P N = ( v (i) - v (2)) p o: 



P f (P) = P^(P) ( 1 - ^ P 2 ) + eoX^ + ^ P 3 + ■ • • . (6) 



The function P/v(P) has a meaning of the extrinsic polarization response of the polydomain 
film to the electric field; i.e. the response, which is produced by the change in volume 
fractions of the two domain states due to the domain wall motions. 

The response of the net spontaneous polarization to the electric field P/v(P) is controlled 
by nonlocal macroscopic phenomena, such as interaction of the domain walls with crystal 
lattice defects [10, 11], depolarizing fields due to the presence of electrode-adjacent layers 
[12, 13, 14], reduction of the domain wall mobility due to free charges [15, 16], etc. In 
order to calculate the function Pjy(P) we adopt the thermodynamic potential G w (Pn, P) 



(per unit area of the ferroelectric capacitor) of the system with the constant electric field 
E. The function G w (Pn, E) includes two contributions: First, the free energy, which is 
increased when the domain wall is shifted from its equilibrium position (e.g. energy of 
the depolarizing field, surface energy of the domain walls, etc.), and, the subtracted work 
produced by the external sources that keeps the constant voltage U = Eh on the electrodes 
(see Fig. 1). In the case of the non-poled (or depolarized) polydomain ferroelectric sample, 
the three lowest terms in the expansion of the function G w (Pn, E) with respect to the net 
spontaneous polarization P N equal: 



G W (P, E) = h 



-ctw-Pjv + ^PwPn ~ PnE 



(7) 



In the case of the stable domain pattern, when a w > 0, the extrinsic polarization response 
is given by the condition for the minimum of the function G w (Pn, E): 

dG w (P N , E)/dP N = (8) 

and the net spontaneous polarization Pn can be expressed in the form of the Taylor series 
with respect to the electric field E: 

P N {E) =e oXw E + ^e b w E 3 + ■■■ , (9) 

where Xw — l/(co^) is the extrinsic contribution to dielectric susceptibility of the system 
and b w = —3/3 w / (eqoc^,) is the nonlinearity constant of the extrinsic contribution to dielectric 
susceptibility of the system with respect to the electric field. 

Substituting Eq. (9) into Eq. (6) one gets the average nonlinear polarization response of 
the polydomain ferroelectric in the form: 

P f (E) = e (Xc + Xw)E+-e (b c + b w ) E 3 + ■ • • , (10) 

where 

c Pi V 8 X J ( j 

is the nonlinearity constant due to the nonlinear dielectric response of the crystal lattice. It 
should be noted that the value of b c is also controlled by the linear extrinsic contribution to 
permittivity Xw Finally, the effective dielectric susceptibility of the polydomain ferroelectric 
layer X f — (^-/ £ o)(9Pf/9E) is equal to: 

Xf(E)^(Xc + X w ) + (bc + b w ) E\ (12) 
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III. METHOD TO STUDY AGEING OF THE POLYDOMAIN FERROELECTRIC 



As it has already been mentioned, we consider that the extrinsic contribution to the 
nonlinear dielectric response \ w + b w E 2 is controlled by a particular nonlocal phenomenon 
associated with some microstructural mechanisms affecting the quality of the ferroelectric. 
It is reasonable to consider that in the systems with the pronounced ageing of ferroelectric 
properties, the associated microstructural mechanism should be also responsible for the 
evolution of the parameters \ w and b w . Therefore, it is reasonable to expect that there 
exists a specific relation between the parameters Xw and b w during ageing, which can be 
expressed in a form Xw — F(bw)- The crucial point of the presented method is the fact, 
that the explicit form of the function F is unique for the every particular microstructural 
mechanism of the dielectric ageing. Hence, when the experimentally observed dependence 
of F(b w ) vs. Xw matches the theoretical prediction, there exists a reasonable indication that 
the considered mechanism is responsible for the observed ageing of the dielectric response. 
It means that this method can be used for the identification of the origin of the ageing of 
ferroelectric films. In this section, we show a general way to determine the characteristic 
function F. 

At first, we shall start with the fact that in the most of the ferroelectric materials with 
pronounced ageing of the dielectric properties, the dielectric nonlinearity is dominated by 
the contribution of the domain wall motion, i.e. b c <^ib w . Then we consider that the evolu- 
tion of the dielectric response is induced by the evolution in some "internal" microstructural 
parameter (e.g. average density of the pinning centers on the domain wall, average do- 
main wall spacing, thickness of the electrode- adjacent dead layer, spontaneous polarization 
screening by free charges, etc.), which we denote by a general symbol £. We consider that 
this parameter controls the values of Xw and b w and we can write: 

Xw = f(0, (13a) 
b w = g(£), (13b) 

and hence 

Xf(E)tt X c+ f "(0+0(0 E 2 . (14) 
If the function g{£) is monotonic, we can express the parameter £ from Eq. (13b), i.e. 
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£ = g 1 (b w ) and substitute it into Eq. (13a): 

X W = f[g- l ib w )] = F(b w ). (15) 

Now, if we consider that the intrinsic contribution to the dielectric susceptibility \ c does not 
change in time, the function F can be cross-checked by the analysis of the data extracted 
from the dielectric ageing measurements on the polydomain ferroelectric films. The first 
measured parameter is the small signal permittivity Sl' 

e L = 1 + X f{0) = £c + Xw, (16) 

where e c = 1 + Xc is the relative permittivity of the crystal lattice. The second measured pa- 
rameter is the dielectric nonlinearity constant b, which can be determined from the following 
expression: 

bE 2 = e f (E)-e L , (17) 

where £f(E) is the electric field dependence of the permittivity of the polydomain ferroelec- 
tric film. 

To apply our method on the experimental data, one should take two steps: first, to 
prove the quadratic dependence of the permittivity on the applied electric field given by Eq. 
(12) and, second, to identify the relationship between the small signal extrinsic permittivity 
Xw and the dielectric nonlinearity parameter b given by Eq. (17) during the dielectric 
ageing. Actually, verifying the first step represents the experimental evidence for a fast 
reversible movement of the domain walls, which should be distinguished from the Rayleigh- 
type relation, where the dielectric permittivity increases linearly with the ac-field as a result 
of the irreversible movement of domain walls under the sub-switching fields (usually a few 
tens of kV/cm) [8, 17, 18]. In the second step, we can identify the particular ageing scenario 
from the experimental evidence of the mutual relationship between the parameters Xw and 
b w in time during the dielectric ageing. By combining Eqs. (15) and (16) one gets 

F{b) = Xw = e L - e c . (18) 

Therefore the validity of Eq. (18) can be demonstrated by the linear relationship between 
the values of F(b) and the total small field dielectric permittivity e^. The relationship given 
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in Eq. (18) can be tested experimentally by fitting the values of F(b) and El to the function 
of the form: 

F(b) = Ke L -B. (19) 

If the considered microstructural mechanism is responsible for the evolution of the dielectric 
response, the fitted value of the parameter K should be approximately equal to unity, i.e. 
K 1, and the fraction B/K should correspond to the value of crystal lattice permittivity, 
i.e. e c ~ B/K, which can be checked by the measurements on the polarized single-domain 
sample. 

IV. EVOLUTION OF THE NONLINEAR DIELECTRIC RESPONSE DURING 
AGEING 

In this section, we demonstrate the applicability of the general method presented above 
for two particular systems. At first we present the system, where the extrinsic permittivity 
is controlled by bending movements of the pinned domain walls. Then we demonstrate 
the system with electrode-adjacent passive layers, where the domain wall movements are 
controlled by the depolarizing field in the electrode-adjacent passive layer. 

A. Ferroelectric ageing in the system with pinned domain walls 

In this subsection, we adopt a model presented in [11], where the extrinsic permittivity of 
the polydomain ferroelectric is controlled by bending movements of the pinned 180° domain 
walls. This scenario seems to be quite reasonable due to several direct observations of such 
interactions of the domain walls with bulk crystal lattice defects [19]. In addition, recent ab- 
initio calculations [20, 21] have indicated that the oxygen vacancy, which is a very common 
defect in the perovskite ferroelectrics, has a smaller formation energy in the 180° domain 
wall than in the bulk and, thus, the oxygen vacancy accompanied with two free electrons 
to comply the requirement of the electroneutrality of the system can be responsible for the 
pinning effect on 180° domain walls. 

Figure 2 shows the model of bending movements of pinned 180° ferroelectric domain wall 
with the attached Cartesian coordinate system. We consider that the pinning centers on the 
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FIG. 2: Model of the 180° domain wall bending. Planar position of the domain wall (doted 
line) is locally blocked by the pinning centers (gray circles). Orientation of vector of spontaneous 
polarization is indicated by gray arrows. When the external electric field E is applied to the 
ferroelectric film, the domain wall is bent between the pinning centers and the profile of the wall 
deflection is described by the function u(x,y). 

domain wall are distributed in the x and y directions with an average distance r between 
them and, thus, the pinning centers density on the domain wall is proportional to 1/r 2 . 
When the electric field E is applied to the ferroelectric film along the ferroelectric axis, the 
originally planar domain wall is bent between pinning centers, as indicated in Fig. 2, and 
the profile of the wall deflection is described by a function u(x, y). The domain wall spacing 
is denoted by a symbol a. 

When the domain wall is bent, the thermodynamic function G w consists of two contribu- 
tions: first, the surface energy of the bent domain wall and, second, the electrostatic energy 
associated with the depolarizing field due to discontinuous change of the normal component 
of the spontaneous polarization at the domain wall. At the constant electric field E, the 
thermodynamic potential G w is given by the formula: 



where the integrals are taken over the square A where x and y are running from —r/2 to 
r/2, a symbol a w stands for the surface energy density associated with the surface tension 
of the bent domain wall, and a symbol (fd stands for the electrostatic potential on the 
domain wall associated with the depolarizing field, which is produced by the bound charges 
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of the surface density due to the discontinuous change of the normal component of the 
spontaneous polarization at the domain wall. The functions u x (x,y) and u y (x,y) are the 
partial derivatives of the function u(x, y) with respect to x and y, respectively. The function 
(fd can be readily calculated by solving the Poisson equation for the electrostatic potential 
and by considering the conditions for the continuity of the normal component of the electric 
displacement and for the continuity of the electrostatic potential. If one considers that the 
net spontaneous polarization P N is equal to: 

2P 



N 



ar 2 



u(x,y)dA, (21) 



and that the maximum deflection u max of the domain wall is much smaller than the average 
distance between the pinning centers r, our detailed calculations [22] show that the leading 
terms in the Taylor expansion of the thermodynamic function G w in this particular system 
with respect to the net spontaneous polarization are: 



G w (Pn, E)^h 



&&w j-,2 0.17 G d w 4 



(22) 



where the symbol a w stands for the thickness of the domain wall. Solving Eq. (8) for this 
particular system, one gets: 

Xw(r, a) w ^f- (23a) 
a a w 

, , . 518.4 r 4 5g £ 3 
a a 3 w Pj 

Now in this particular system, we can distinguish two scenarios of the dielectric ageing. 



1. Dielectric ageing due to the progressive pinning 

In this scenario, we consider, that due to the drift of the crystal lattice defects (e.g. oxygen 
vacancies) in the material, the concentration of the pinning centers on the domain wall may 
be increasing during the dielectric ageing. Therefore, in this scenario, which is characterized 
by a progressive pinning of the domain walls, we consider that the microstructural parameter 
responsible for the dielectric ageing is the average distance r of the pinning centers and that 
the domain wall spacing a remains constant. 
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Then we can express the parameter r from Eq. (23b) and by substituting it into Eq. (23a) 
we arrive at the following form of the function F(b): 



F(») = VS.teS(2=y (24) 



In addition, from the time evolution of the parameters b(t) and we can estimate the 
evolution of the average distance r(t) between pinning centers on the domain wall: 



= 0.088 a „P b(t) 

2. Dielectric ageing due to the domain wall coalescence 

On the other hand, the same microstructural model offers a different ageing scenario. In 
the ferroelectric polydomain film, it is possible that the domain pattern is in the essentially 
non-equilibrium configuration and, therefore, there exists a "thermodynamic" force applied 
to the domain wall that drives the system to reach the "absolute" equilibrium configuration. 
In this situation, the depinning of some walls can occur and they can coalesce to reduce the 
energy proportional to the domain wall area. Then, we can consider that the parameter 
responsible for the evolution of the dielectric response is the domain spacing a and that the 
concentration of the pinning centers on the domain wall remains constant during ageing. 
Hence, we can express the parameter a from Eq. (23b) and by substituting it into Eq. (23a) 
we obtain the characteristic function F(b) in the form: 

2 



F(b) 



0.088 P n fa 



(26) 



e e c \ r 

Similarly as in the previous ageing scenario, it is possible to estimate the time evolution of 
the average domain spacing a(t) from the time evolution of the parameters b(t) and e^t): 

= 0.031 a.FSHt) 
e\e c [s L (t) - e c ] 

Thus, it is seen that using cross-checking the two qualitatively different b vs. El depen- 
dences, it is possible to distinguish the microstructural processes, which are responsible for 
the evolution of the nonlinear dielectric response. 
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FIG. 3: Scheme of the electroded ferroelectric film with passive layers. It is always considered 
that h » d. 

B. Ferroelectric ageing in the system with the electrode-adjacent passive layer 

In this subsection, we study the dielectric response of the system, where the ferroelectric 
layer with the antiparallel domain pattern is separated from the electrodes by the electrode- 
adjacent passive layers. Although the nature of these layers is not clear [23] and the direct 
observations of well-defined layers are rather scarce [24], numerous indirect experimental 
observations support the wide acceptance of the passive layer concept. Until today, various 
mechanisms of film-electrode interactions have been proposed: creation of variable depletion 
layers associated with the Schottky barriers [25, 26], low permittivity layers due to the 
oxygen vacancy entrapment at ferroelectric-electrode interfaces [27], local suppression of 
polarization states in the near-electrode regions [28] , local diffusion of the electrode material 
into the ferroelectric [29], or a chemically distinct surface phase [30]. 

In order to calculate the dielectric response of the system with a passive layer, we adopt 
the model shown in Fig. 3. We consider a sandwich structure of the ferroelectric layer of 
the thickness h with a lamellar antiparallel domain pattern, two dielectric passive layers of 
the total thickness d and the permittivity Ed, which separate the ferroelectric layer from 
the electrodes. The spontaneous polarization Pq is considered to be constant within each 
domain and its abrupt change to zero at the interface of the ferroelectric and dielectric layers 
[z = ±/i/2) leads to the appearance of a bound charge of the surface density at, = Pq. Due 
to finite conductivity of the passive layer, the bound charge is considered to be compensated 
by the free charge injected from the electrodes, so the proportionality factor A called the 
degree of screening is introduced to measure the process of the charge transport er/ = — \o~b, 
where < A < 1. 
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In this model, we neglect the effect of the lattice pinning of the domain wall and we 
consider, for simplicity, that the net spontaneous polarization response to the electric field 
is controlled by the electrostatic energy of the depolarizing field produced by the periodic 
distribution of the bound o"b and free <7/ charges at the interface of the ferroelectric and 
dead layers. In addition, we address in detail the case of the neutral domain pattern, i.e. 
we set P/v = in the absence of the applied voltage U = 0. The thermodynamic potential 
G w (per unit area of the sample) of a system with a given voltage U on the electrodes and 
free charges <jf inside the system, which includes the electrostatic energy of the electric field 
with the subtracted work produced by the external voltage sources, is given by the formula: 

G w = / Ei (A - Po,i) dV + Ua , (28) 



2Ajy 

where A is the area of the ferroelectric capacitor and <To is the surface density of free charges 
on the electrodes. 

Calculation of the electrostatic energy requires the knowledge of the electrostatic potential 
tp, which satisfies the following equations of electrostatics in the ferroelectric and the passive 
layers, 

£a ^Pi + £c ^Pi = Q ^ #Vp + ^V, =0; (29) 

a dx 2 c dz 2 ' dx 2 dz 2 
with the boundary conditions <p = —(+)U/2 at z = +(—)(h + d)/2, and 

£ ^- £ ^ = ^p <Pf = <p d , atz = /i/2, (30) 

where the subscripts / and g denote the ferroelectric and passive layer, respectively. Here 
a(x) = crfe(x) + (Tf(x) is the sum of the surface density of the bound charge due to the 
spontaneous polarization and the free screening charge from the electrode. Solution of the 
potential and the calculation of the electrostatic energy for a similar model are already 
available in the literature [12, 13, 14], however with a different distribution of charges at 
the interface of the ferroelectric and passive layers. In this model, we consider that the 
spontaneous polarization continuously and smoothly changes its value within a distance a w 
representing the domain wall thickness. Then the function a(x) can be obtained in a form 
of the Fourier series: 



^ 64a 2 P . 2 /nira w \ ( 

n=l 



nixx 



n-K ( P N 
2 V Po 



X sin ( 
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(31) 



In order to calculate the dielectric response, we express the function G w in terms of the 
applied electrode voltage U and the net spontaneous polarization of the ferroelectric layer 
Pn- With use of the results published earlier [12, 13, 14], the leading terms of the function 
G w in this model are equal to: 



G w (Pn, E) 

where 



- ®wPn + t PwPn ~ PnE 



and 

feln-irk 



niiT / [e^i 

D n {T) = e d coth + V^lZc coth I J — - 



(32) 



f\ A - _ I 1 - 67fl5 A ~ (~ 1 )"( A ~ 2 ) • 4 ( ™a w \ 

aw{ ' J ~ e (e d h + e c d) + e^Ji ^ n*D n (d/a) Sm I 4a i' ^ 



0.688 a 5 ^ (-l)"(A-8)-A . 4 /nvra^x 
^ (A ' rf) = g n3D n (^/a) Sm (^-) ' (33b) 

J5 = f dU (33c) 



(33d) 



2a 

Similarly as in the previous section also in this particular system we can distinguish two 
ageing scenarios: 



1. Dielectric ageing due to the thickening of the electrode- adjacent passive layer 

It was reported [27, 31] that the fatigue of the dielectric response of the ferroelectric films 
can be caused by an increase of the thickness of the passive layer (e.g. due to the oxygen 
vacancy entrapment at the film-electrode interface). The main indication, from which such a 
statement is inferred, was the measurement of the size effect on the permittivity. There there 
was observed a decrease of the effective capacitance of the interface layer. Unfortunately, it 
was also demonstrated [15] that the same effect on the permittivity, i.e. the effective decrease 
of the interface layer capacitance, can be produced by the spontaneous polarization screening 
by immobile charged defects. Nevertheless, by applying the present method on the data of 
the nonlinear permittivity, it was possible to distinguish whether such a fatigue scenario 
may take place. 
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In order to make it possible to use the principles of the discussed method to this particular 
situation, we should make following considerations: (i) The dominant phenomenon, which 
controls the extrinsic dielectric response, is the increase of the passive layer thickness, d > 0; 
(ii) There is no spontaneous polarization screening, A = 0; (iii) The domain pattern in the 
ferroelectric layer is dense, a -C h; (iv) The passive layer in the system is thin compared 
to the average domain spacing, d <C a; (v) The thickness of the domain wall is much 
smaller than the domain spacing, a w <C a; (vi) The permittivity of the passive layer can 
be approximated by the permittivity of the ferroelectric in the following way: = ^s a e c ; 
Finally (vii), the average domain spacing a remains constant during the dielectric ageing. 
Under these considerations, Eqs. (33) can be further simplified and Eq. (32) has a simple 
form: 



G W (P N , E)^h 



F N + 1QOn . .p2 /= F N - F N F 



(34) 



_8ahe 0y /e^r c 192 ahe P^ \yfe^e c 
Solving Eq. (8) for this particular situation, one gets: 

Xw{d) « , (35a) 

dr 7r 

^ 16a 3 fe 3 eg {e a e c f 2 
b w (d) ~ . (35b) 

Finally, the function F(b) is obtained by expressing the parameter d from Eq. (35b) and by 
substituting it into Eq. (35a): 



F{b) = ^bf-^l. (36) 



s 



o 



2. Dielectric ageing due to the spontaneous polarization screening 

Here we show that within the same microstructural model of the passive layer, we can 
discuss the possibility of a different ageing scenario. Due to the essentially semiconductive 
nature of ferroelectric materials and due to the strong experimental evidence [32, 33], the free 
charge injection across the electrode- adjacent layer has become a widely accepted concept. 
Nevertheless, it has been discussed mainly in connection with the ferroelectric switching 
measurements and the phenomenon called imprint, i.e. the horizontal shift of the hysteresis 
loop, which is caused by the internal bias field produced by the injected free charges. In this 
subsection, we deal with the free charge injection across the electrode-adjacent layer in the 
system with the polydomain ferroelectric layer. 
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It has been already reported [31] that the pinning of domain walls in the electrode- 
adjacent region may represent a source of the essential nonlinearity of the dielectric response. 
It has been also reported [15] that the similar pinning effect on the domain walls can be 
produced by the low mobility of the free charge at the film-electrode interface. Here we 
study how the injected free charges affect the nonlinear dielectric response of polydomain 
films. 

In order to make a clear statement, whether the free charge injection across the passive 
layer can be responsible for the observed dielectric ageing using the presented method, we 
shall make the following considerations: (i) The dominant phenomenon, which controls the 
extrinsic dielectric response is the spontaneous polarization screening, A > 0; (ii) The domain 
pattern in the ferroelectric layer is dense, a <C h; (iii) The passive layer in the system is thin 
compared to the average domain spacing, d <^ a; (iv) The thickness of the domain wall is 
much smaller than the domain spacing, a w <C a; (v) The permittivity of the passive layer 
can be approximated by the permittivity of the ferroelectric in a following way, Ed = y/£ a Ec] 
(vi) Since we are focused on the deteriorating impact of the free screening charges on the 
dielectric response of the polydomain ferroelectric film, which is produced due to the pinning 
effect of the free charges on the domain walls, we further consider that the free charges are of 
much lower mobility than the mobility of the domain wall in the absence of the free charges; 
Finally (vii), the average domain spacing a remains constant during the dielectric ageing. 
Under these considerations, Eqs. (33) can be further simplified and Eq. (32) has a simple 
form: 



G w (Pn, E)^h 



2 A ad 2 A a d 



- ~ c ud2 r^ = ~~ PnE 

_3e a w hy/e a e c §e Q a%hP$y/e a e c 
where E = U/h. Solving Eq. (8) for this particular situation, one gets: 



(37) 



/x x 3 a w h y/ E a e c 

XUX) 4adX ' (38a) 

M A) „ °- 633 ;f A ^ )3/2 . (38b) 

It should be stressed again that the limit for A approaching zero cannot be applied to 
Eqs. (38), since it violates the conditions for the applicability of our approximations. The 
reason for it can be described as follows. When the degree of screening A is smaller than some 
critical value A cr ;t, the leading terms that control the thermodynamic function G w given by 
Eqs. (32) and (33) become very weakly controlled by the spontaneous polarization screening 
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A and the other parameters of the system have stronger effect on the net spontaneous 
response of the system. The critical value A cr it can be easily calculated by comparing the 
quadratic term coefficients in the expansion of the thermodynamic function G w with respect 
to the net spontaneous polarization in the two aforementioned different approximations given 
by Eqs. (34) and (37); this value reads: A cr i t ~ a w d/a 2 . 

When we express the parameter A from Eq. (38b) and when we substitute it into Eq. (38a), 
we obtain the characteristic function F{b) in the form: 



In addition, it is possible to estimate the time evolution of the process of charge injection 
measured by the degree of screening A(t) from the time evolution of the parameters b(t) and 



However, it is seen that the thickness of the passive layer d enters into the above equation 
and its value is always difficult to determine experimentally. Instead of that, we can use the 
expression for the effective capacitance per unite area of the ferroelectric-electrode interface 
Cd ~ toy/Za^c/d, which can be estimated more easily. 

V. DISCUSSION AND CONCLUSIONS 

We have theoretically analyzed the nonlinear dielectric response of the ferroelectric poly- 
domain film in two systems with different microstructural features: first, the system with 
the pinned domain walls and, second, the system with the elect rode- adjacent passive layer. 
In these two considered systems, we studied four possible scenarios of the dielectric aging 
and found that there exists a clear qualitative difference in the evolution of the nonlinear 
dielectric response in each particular scenario. We found that the relation between the 
nonlinearity constant b w and the small signal permittivity El expressed by the function 
F(b w ) = Xw — £l — ^ci which is shown in Tab. I, is different in the each particular ageing 
scenario. Therefore, the analysis of the experimentally observed evolution of the dielectric 
nonlinearity constant b and the small signal permittivity El can be used as a method to 




(39) 



e L (t): 




(40) 
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identify the micro-structural mechanism, which plays the dominant role in the evolution of 
the nonlinear dielectric response. 

To identify whether the particular microstructural mechanism is dominantly responsible 
for the dielectric ageing, the experimentally measured values of b and El should be fitted to 
the function of the general form: 

F(b) = Ke L -B. (41) 

If the considered mechanism controls the ageing, the fitted constant K should approximately 
equal to unity and the constant B should be equal to e c . 

However, applicability of the presented method is subject to several considerations that 
should be kept in mind: (i) We consider that the intrinsic permittivity e c is time indepen- 
dent, (ii) We limit our analysis to the systems of uniaxial ferroelectric film with the 180° 
domain pattern and we neglect the effect of the non-180 domain walls on the dielectric 
response, (iii) We consider that the domain pattern in the ferroelectric is neutral, it means 
that the net spontaneous polarization is zero in the absence of the electric field, (iv) We 
consider that the nonlinear dielectric response is dominated by the domain wall movement 
and not by the dielectric response of the crystal lattice, b c -C b w . The first assumption is 



TABLE I: Characteristic functions F(b), which relates the constants of dielectric nonlinearity b 
with the small signal extrinsic contribution to permittivity Xw of the polydomain ferroelectric films. 



Ageing scenario 



F(b) 



Progressive pinning by crystal lattice defects Vb / p o 



Domain wall coalescence 



0.088 Pp (q_ 

e e 



Thickening of the electrode-adjacent layer 



0.41 P? 



Spontaneous polarization screening 
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intuitive. To justify the second one, we need to point out that the extrinsic contribution to 
the nonlinear permittivity of the non-180 domain walls is negligible due to the elastic effects 
suppressing their mobility. The third assumption puts the restraints on the preparation of 
the ferroelectric samples, which should be depoled, e.g. using a fast decayed low frequency 
ac field. Finally, the last assumption should be checked for each ageing scenario and the 
particular configuration of the domain pattern microstructure using the fitted experimental 
values. 

We believe that the presented method can by used as a simple and useful tool for getting a 
deeper insight into the configuration of the domain pattern and the quality of the ferroelectric 
thin films and for providing a way to identify the dominant microstructural phenomenon, 
which is responsible for the ageing of dielectric response of ferroelectric polydomain film. 
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